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We theoretically study an extension of the Dicke model, where the single-particle Hamiltonian has 
three energy levels in Lambda-configuration, i.e. the excited state is coupled to two non-degenerate 
ground states via two independent quantized light fields. The corresponding many-body Hamiltonian 
can be diagonalized in the thermodynamic limit with the help of a generalized Holstein-Primakoff 
transformation. Analyzing the ground-state energy and the excitation energies, we identify one 
normal and two superradiant phases, separated by phase transitions of both first and second order. 
A phase with both superradiant states coexisting is not stable. In addition, in the limit of two 
degenerate ground states a dark state emerges, which seems to be analogous to the dark state 
appearing in the well known stimulated Raman adiabatic passage scheme. 
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I. INTRODUCTION 

Superradiance is a collective phenomenon originating 
from atomic physics. There, it is regarded as a collec- 
tive spontaneous emission process of a dense ensemble of 
radiating atoms [1 . The atoms interact indirectly via a 
light field. The first microscopic description of this phe- 
nomenon was given by Dicke |2J. 

In the context of phase transitions a collection of two- 
level systems coupled linearly to one scalar bosonic mode 
undergoes a second-order phase transition from a nor- 
mal to a superradiant phase at a certain critical coupling 
strength. This phase transition has been investigated 
theoretically a long time ago by Hepp and Lieb [3] and 
also by Wang and Hioe However, there is no ex- 
perimental realization in atomic systems to date. There 
were theoretical proposals to produce this phase transi- 
tion in artificial quantum systems like circuit or cavity 
quantum electrodynamic (QED) systems [5-7 . Though, 
there exist no-go theorems for atomic, cavity and circuit 
QED systems which theoretically preclude the normal- 
superradiant phase transition [8-10_. 

Recently, experimental progress was achieved in this 
field by the group of Esslinger, who coupled a Bose- 
Einstein condensate to a single mode of an open opti- 
cal cavity [TT]. The unitary dynamics of this system is 
described by an effective Dicke Hamiltonian |111I12) . Ex- 
perimentally, the normal-superradiant phase transition 
is observed by measuring the mean intracavity photon 
number. 

Inspired by this experimental realization of an effective 
Dicke-Hamiltonian, in this paper we theoretically inves- 
tigate an extension of the Dicke model. Here, three-level 
systems in A-configuration are considered, which are cou- 
pled to two independent scalar bosonic modes. We are 
interested in how the phase transition is changed in this 
configuration. Furthermore, coherent population trap- 
ping [13], dark states and the stimulated Raman adiabatic 
passage (STIRAP) scheme [14J are associated with this 
kind of system in the single-particle and semi-classical 
case. We therefore study to what extent dark state 




Figure 1: (Color online) Level structure of the A- 
configuration. One particle has two ground states |1), 2) 
and one excited state |3), where the excited state is coupled 
to the two ground states via two independent scalar bosonic 
modes with in general different frequencies u>2 and cou- 
pling strengths gi, g-2. 

physics plays a role in our quantum-many-body setting. 

The paper is organized as follows: At the beginning, 
in Sec. [TT] we introduce the model, give a detailed de- 
scription of the Hamiltonian and discuss the symmetries 
of the model. Subsequently in Sec. |III| we describe the 
Holstein-Primakoff transformation for multilevel systems 
and derive an effective Hamiltonian in the thermody- 
namic limit. We diagonalize this effective Hamiltonian 
and give explicit expressions for the ground-state energy 
and the excitation energies. Section |IV| addresses the 
phase transition: The zero-temperature phase diagram 
is mapped out and analysed. In Sec |IV B| we discuss 
properties of the appearing dark state. Finally, Sec. [V] 
closes our contribution with some conclusions. 



II. THE MODEL 

We consider a quantum mechanical system consisting 
of M distinguishable particles and two independent scalar 
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bosonic modes. Each particle i possess three energy levels 
|l) (i) , |2) (i) , and |3) (i) with energies E 1 < E 2 < E 3 , 
respectively. For later analysis we define, A = E3 — Ei, 
6 = E 2 - Ei, with A > 5 > 0. The level scheme is 
in so-called Lam&da(A)-configuration. Each of the two 
lowest energy levels couple to the highest energy level 
via one of the bosonic modes respectively (see Fig. [T]). 
The Hamiltonian has the form (h = 1) 
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H = ^E n A n n + {^n &l a n 
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+ ^(a n + a n )(A*+Af)}. (1) 

Here, A r s are defined by 

AT 

V = 5>> W < S I' r, S = l,2,3, (2) 

i-l 

and represent collective particle operators. 

The diagonal operators A™ measure the occupation of 
the nth energy level, i.e. how many of the M particles 
are in the energy state \n). This illustrates the first term 
in the Hamiltonian 0. The second term gives the en- 
ergy of the two scalar bosonic modes, each one having 
the excitation energy u>i and ui2 respectively. The op- 
erators a n and a n create and annihilate a boson in the 
nth mode. They fulfill canonical commutator relations, 
[a n ,aln\ — 6 n>m and [a n ,a m ] = 0. Lastly, the third term 
in the Hamiltonian represents the interaction of the 
particles with the two bosonic modes. Here, g n are the 
corresponding coupling constants. 

We call the first |1) and the third |3) energy levels of 
the particle system together with the first bosonic mode 
blue branch, since A > 5 is assumed. Correspondingly, 
we call the second |2) and the third |3) energy level of the 
particle system together with the second bosonic mode 
red branch (cf. Fig. [TJ . 

Symmetries and phase transition. Our Model is a gen- 
eralization of the Dicke model [2[T5], where particles with 
only two energy levels are considered, and the two states 
are coupled via one scalar bosonic mode. In the thermo- 
dynamic limit the Dicke model exhibits a non-analytic 
behavior in physical observables as a function of the cou- 
pling strength g. Thus, the Dicke model exhibits a (quan- 
tum) phase transition, which is continuous, i.e. of second 
order and separates two phases: a normal phase and a 
so-called superradiant phase. The superradiant phase has 
a ground state with spontaneously broken symmetry. A 
similar behavior is anticipated in the extended model. 

In analogy to the Dicke model, here exist two symme- 
try operators 

n n = exp\-iTt(-A n n + a n a n )\, n=l,2, (3) 



which commute with the Hamiltonian ([T]). These oper- 
ators have the physical meaning of parity operators and 
have eigenvalues ±1. The operator fj n = —A n n + a n a n 
in the exponent of the parity operator ^ is related to 
the number of excitations in the blue (n — 1) or in the 
red (n = 2) branch of the A-system and the number of 
excitations in the corresponding nth bosonic mode, re- 
spectively. The operator fj n itself is not conserved, i.e. 
[fj n , H] ^ 0. This is consistent with the Dicke model |15| . 
In the rotating wave approximation the operators f) n be- 
come conserved quantities. Conservation of the two par- 
ities means that the Hilbert space decomposes into four 
irreducible subspaces. It is the parity which is sponta- 
neously broken in the superradiant phase of the Dicke 
model. Thus, we expect that at least one of the parities 
is also spontaneously broken in our model. 

Using the definition ^ of the operators A r s , one 
can show that the two sets of operators {|(^4.3 3 — 

A^), An, ^3™}; n = 1,2 fulfill the angular momentum 
algebra respectively, i.e. they are generators of the spe- 
cial unitary group SU(2) and can be understood as an- 
gular momentum operators. In addition, the operators 
A r s fulfill the algebra of generators of the unitary group 

u(3) nnun] 

[A r , A n ] = 5 s ,n A r — <5 r , m A n (4) 

and are, according to that, generators of the group U(3). 
It is known, that the generators of the group U(N) can 
be represented by either N or by N — 1 independent 
bosons |16l [T7] . The first choice corresponds to the 
Schwinger boson representation |16] , the latter choice 
to the Holstein-Primakoff transformation of the genera- 
tors [MIH]- 

III. METHODS 

The Dicke model was introduced in 1954 [2 - To date 
there exists no exact analytical solution to this model 
for a finite number M of particles. However, the Dicke 
Hamiltonian can be exactly diagonalized in the thermo- 
dynamic limit [15 , i.e. Af — > 00. This can be achieved by 
using the already mentioned Holstein-Primakoff trans- 
formation. In this article we apply a generalized version 
of the Holstein-Primakoff transformation to diagonalize 
the Hamiltonian ([lj of the A-system. 

A. Generalized Holstein-Primakoff transformation 

In the present paper we discuss the A-system which has 
N = 3 single-particle states. Though, we will formulate 
the following argument for a general number N of single- 
particle states. The number of particles is denoted by Af, 
whereas the number of single-particle states is denoted by 
N. 



The generalized Holstein-Primakoff transformation 
maps the generators A r s of the group U(N) onto a com- 
bination of creation and annihilation operators b r ,b r of 
N — 1 independent bosons. Hence, the operators b\. and 
b r fulfill canonical commutator relations, [b,.,bl] = 5 rs , 
[b r ,b s ] = 0. These bosons we will refer to as Holstein- 
Primakoff bosons (HP bosons). One of the N states of 
the one particle system is called the reference state, which 
we denote with |m). The meaning of the state |m) and 
which of the N states can be used as a reference state 
will be elucidated later. Then, the generalized Holstein- 
Primakoff transformation is given by [17] 
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Figure 2: (Color online) Physical interpretation of the bosons 
introduced via the generalized Holstein-Primakoff transfor- 
mation (JHJ: The two bosonic operators bj., with r 7^ m, can 
be understood as collectively exciting the particles from the 
reference state \m) {left: m = 1, right: m — 2) to the state 
r). The analogue holds for the annihilation operators b r . 



A r s = btb s , 



i r m = Ste m (7V), 
A^ = e m (Af)b s , 

A™ = G m (N) 2 , 



> r, s =/= to 



(5) 



with 



(6) 



There are at most Af HP bosons per mode, i.e. the ex- 
pectation value satisfies (b\ b r ) < Af, r 7^ to, due to the 
operator Q m (Af) and the fact that b r acting on a state 
with zero HP bosons in the rth mode equals to zero. In 
addition, the number of HP bosons in all N — 1 modes 
does not exceed Af, i.e. J2 r ^m V>r M < A/". 

We now apply the generalized Holstein-Primakoff 
transformation ^ to the Hamiltonian (JT|) with e.g. |1) 
as the reference state (to = 1) and obtain 



H m =i = E 1 Af + S b\ b 2 + A b\ 63 + u n a„ a v 



+ 



9i /*t 



aj+ax) (ble 1 (Af)+e 1 (Af)b 3 ) 

— + (bth + bih). (7) 



+ 



Af 



The first line is the free part of the Hamiltonian, from 
which one can infer the meaning of the HP bosons: The 
number of HP bosons in the mode with frequency 5 is 
given by the operator b\ 62 . This means that b\ is related 
to the creation of excitations with energy 5, which is the 
energy separation of the single-particle energy levels |1) 
and 1 2). Thus, the operator b\ can be understood as 
collectively exciting the particles from the first energy 
level to the second one. This is visualized in Fig. [2] An 
analogous reasoning can be given for the other HP boson 
corresponding to the operator 63. 



B. Thermodynamic limit 

The expectation value of the HP boson operators b r 
is zero for a finite number Af of particles. In contrast, 
in the thermodynamic limit the expectation value of this 
operator can be finite, and is then macroscopic. Given 
that the occupations (A„ n ), and (ajj a„) should scale with 
the particle number Af, we make the ansatz 



b r = VAf <F r + d r , m, 
a n = VAf (p n + c n , n 



1,2, 



(8a) 
(8b) 



in the thermodynamic limit. Here vA^ ^ r an d VAf tp n 
are the ground-state expectation values of b r and a„, re- 
spectively. This means that the ground-state expectation 
value of the bosonic operators d r and c„ is zero and, con- 
sequently, these operators can be interpreted as quantum 
fluctuations. Furthermore, they fulfill canonical commu- 
tator relations and their matrix elements are of the order 
of Af°. The parameters ^ r and (p n can be chosen real 
and range from zero to one, which ensures {b\b r ) < Af. 
Another viewpoint is, that the operators d r and c„ can 
be generated from b r and d„ respectively by a canon- 
ical transformation and can be considered as displaced 
bosonic modes |15j . 

Using the ansatz Q we find that the ground-state oc- 
cupations of the particles and of the scalar bosonic modes 
are given by 



(A r r ) =Afy 2 r + (did r ) , r^m, 
(A™)=N^ 2 m -Y(did r ), 



(a+ On) =Aftfl+ (4c«) 

with the abbreviation 



n=l,2, 



V m = 1 



(9) 
(10) 

(11) 



(12) 



Inserting the ansatz Q into the operator O m (A/ - ) ^ 
of the Holstein-Primakoff transformation ^ we obtain 
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C. Ground-state properties 



9 m (W) = Nlprn - £ + * r (4 + d r )] . 

. (13) 

Since we are working in the thermodynamic limit, we 
can asymptotically expand the square root in powers of 
Wl/Af and obtain up to the order AT -1 : 



e m (A0 « vOv Vv, 
i 



1 



£*,(4 



2A/> 2 



£[44+£^(4+4)(4+4) 



(14) 

In this expansion we have neglected terms of the order 
7V _3//2 and higher, which do not contribute to the Hamil- 
tonian (7]) in the thermodynamic limit. 

Finally, we insert the expression ( 14 1 for the operator 



m (A/") and the ansatz ^ into the Hamiltonian In 
the thermodynamic limit we can neglect terms with in- 
verse powers of Af and constants of the order Af° . This 
eventually yields 



h 



(2) 



with 



n m=1 - Ex 



6*i 



L m=l ~ 



(15) 

(16) 
(17) 



■Agnp\ Ipl * 3 + 4 #2 y>2 *2 * 3 , 

+ 2 ff2 <^2 * 3 ] + (4 + 4) [ A *3 

+ 2 5l ip 1 ipi (1 - *3/t/> 2 ) + 2 52 </>> *2] 

+ (cj + ci) (wi tpi + 2 .gi Vi * 3 )] 

+ (4 + C 2 ) (w 2 <y5 2 + 2 32 *2 * 3 )] , 

&£Li=4<fe[*-20m*3M] (18) 
+ ^<fe [A -2ffi ^1*3/^1] 
+ cji c{ ci + w 2 4 c 2 

- (4 + ^2)^31^1*2 *sM 8 

- (4 + 4) V ¥>i * 3 M(i + |*iM 2 ) 

- (4 + 4) (4 + 4)si * 2 m (1 + * 3 M 2 ) 

+ (4^2+44)232 ¥>2 

- (cj + Cl) (4 + 4)31 *2 * 3 M 

+ (4 + ci) (4 + 4)51 ^1(1 - *i/v> 2 ) 

+ (4 + C 2 ) (4 + 4)52 *3 
+ (4+C 2 )(4+4).92 *2- 

The Hamiltonian H separates into three parts h 1 - " ' , each 
one scaling with J\f( 2 ~ n )/ 2 and containing products of n 
operators d r , c^. 



The ground-state energy hm ( |16| of the Hamilto- 
nian (15 1 is a function of the parameters ipx, tf2, *2 and 
^3. Next, we extremize the ground-state energy with 
respect to theses parameters, i.e. we stipulate 



dh {0) . 

un m=l _•_ 

ggli ! 



0, n=l,2, 
0, r = 2,3. 



(19a) 
(19b) 



In the case of ipx being finite, this stipulation is equivalent 



to set the coefficients of the linear Hamiltonian h^Li \Y1 \ 
equal to zero (cf. Ref. |15|V 

The first set of Eqs. ( |19a l gives conditional equations 
for the parameters ip n of the scalar bosonic modes, 



ipx = -2— tpi 



<p 2 = -2— * 2 *3, 
W 2 



(20) 



which, when inserted into the second set of Eqs. (19b I 



gives conditional equations for the parameters ty T of the 
HP-bosons, 



\L1J1 



U>2 



*! 



2 2 
A - 4— fl - *§ - 2*^ - 4^ 2 -1' 2 . 



wi 



= 0, 

*3 =0. 



(21a) 
(21b) 



These equations have several sets of solutions: 

(?) Normal state. The trivial solution, \&2 = * 3 = 0, is 
attended by ipx = tp 2 = (cf. Eq. (20 1). Since </> 2 mea- 
sures the macroscopic (~ Af) ground-state expectation 



value of the nth scalar bosonic mode (cf. Eq. (Ill), this 



trivial solution describes the normal state, i.e. no super- 
radiant state of the system. In addition, the ground-state 
expectation value of the occupation of the nth energy 
level, which is given by (A n n ), is macroscopic for 11 = 1 
only (cf. Eqs. Q, ( 10 1 ) . Thus, all particles occupy their 
respective ground state |1). The ground-state energy of 
the many-particle system is given by h 



(0) 

normal 



E x . Fi- 
nally, we note that the normal state is always a solution 
of the equations (20) and (21 1, irrespective of the cou- 



plings 31 and 32. However, analyzing the Hessian matrix 
of /i^Li restricts the range of the first coupling to 



.91 



< y/Aux/2 = 



.91, 



(22) 



(ii) Blue superradiant state. The second solution of 
Eq. (21 ) is given by 



* 2 = 0, *, = ±a/-W1- 



9i,c 
9i 



(23a) 



_^.9i L 9i,c 
wi V V 31 



<p 2 = 0. (23b) 



5 



In contrast to the previous solution, this solution has a 
finite parameter ipi and for this reason a finite and macro- 
scopic occupation (a\ fii) of the first scalar bosonic mode. 
This solution corresponds to a superradiant state of the 
system, where superradiance occurs in the blue branch 
of the A-system. More precisely, we call this state a blue 
superradiant state. Furthermore, the first and the third 
single-particle energy level are macroscopically occupied. 
If we insert the solution (23 1 into the ground-state en- 



ergy ( 16 1 of the many-particle system, we obtain 



£(0) _ rp A ( 9i 



1 - 



(-y 



(24) 



Hence, the ground-state energy of the superradiant state 
is always smaller than the ground-state energy of the 
normal state. However, this solution is only valid for 
couplings gi > since for smaller couplings gi the 

non-zero parameters of the solution ( 23 1 become purely 



imaginary and, in addition, the Hessian matrix of ft^Li 
becomes indefinite. 

(in) Red superradiant state. There can be another set 
of parameters tp n , ^ r which extremize the ground-state 
energy hm . This set cannot be deduced from the ground- 
state energy h^L~i from Eq. (16 1, because it represents 



not a local but a global minimum of fi^ =1 . Since h^ =1 is 
defined on the unit ball B 2 = {(x, y) € R 2 \ x 2 +y 2 < l}, 
the global minimum lies on the boundary of B 2 , that 
is + = 1 (ipi = 0) holds. To obtain this global 
minimum one has to first of all set '01=0 in Eq. (16) 



and omit all terms involving tpi in Eq. (17 1. Secondly, one 



extremizes the ground-state energy as before, but taking 
the constraint ^ 2 + ^ 2 = l into account. Eventually, we 
obtain 



*2 = ± 




where we have introduced 



92, 



\/(A - S) ui 2 



(25a) 
(25b) 
(25c) 

(26) 



The occupation of the first energy level 1 1) is not macro- 
scopic, i.e. it is negligible in the thermodynamic limit. 
Since <p\ = and ip 2 is finite, this state also corresponds 
to a superradiant state, whereat superradiance occurs in 
the red branch of the A-system. We call this superradiant 
state a red superradiant state. 

This solution can also be found by direct extremiza- 

,(o) " 



considers the second level |2) as the reference state m of 
the Holstein-Primakoff transformation ([5]). In general, 
one can say that using the Holstein-Primakoff transfor- 
mation in the thermodynamic limit with the mth state 
as the reference state, one can describe many-particle 
states in which the occupation of the mth energy level of 
the single-particle system is finite. In order to describe 
the normal state, which is a state where all particles oc- 
cupy their respective ground state |1), one has to take 
|1) as the reference state (m =1). In contrast, to de- 
scribe a state where no particle occupies its respective 
ground state |1), either |2) (m = 2) or |3) (m = 3) 
has to be chosen as the reference state. Furthermore, 
we note that one can easily obtain h^L 2 (^i, ^3) fr° m 



h^' =1 (^i 2 ,^ 3 ) from Eq. (|16j) by substituting ipi i-> *i 

and * 2 ^ ^2 = \/l - *f - *| . 

At last, the ground-state energy of this red superradi- 
ant state is given by 



- TP 

"red - ^1 



-1 2 



92 

92,c 2 



92 



where 



92.. 



A + V5) 



(27) 



(28) 



is a second critical coupling strength. 

Unphysical solution. There is also a solution of the 
Eqs. (211 which corresponds to a state where both 



branches of the A-system are superradiant. However, this 
state is either not well defined for certain couplings g\ and 
g 2 or it does not minimize the ground-state energy (16 1. 
In the latter case, this solution can be attributed to a 
point of inflection on the energy landscape h!^ > =1 (^ 2 , ^3) 



A further solution of the Eqs. (211 represents a dark 



state. This state is discussed in detail in Sec. IIV Bl 



D. Excitation energies 

So far, we have extremize the ground-state energy 
of the Hamiltonian (151 in the thermodynamic limit. By 



this procedure the linear part hP~> is eliminated as well. 
The next step is to diagonalize the quadratic part hP>. 
This can be achieved by means of a principle axis or Bo- 
goliubov transformation |15j . The diagonalized Hamilto- 
nian is then given by 



h {2) = y^£ fc e[ efci 

k 



(29) 



where e\.(ek) create (annihilate) quasi-particles which re- 



tion of the ground-state energy ^=2(^15 ^3)1 l - e - ^ one -^ er ^° bosonic excitations, i.e. 



el and ik satisfy canonical 
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commutator relations. The operators el, ej, and the exci- 
tation energies Ek have to be evaluated separately in the 
three different states. The determination of these quanti- 
ties reduces to a diagonalization of two-by-two matrices. 
The diagonalization procedure yields four excitation en- 
ergies, given by (fc = (x e {1, 2}, a e {+, -})) 



£ x,± — ~ 



2 
± 
1 



uj\ + uj'i _ + 2Awi i 



(30) 



2Awi,_) + 16gju;iu;i 



2 2 
Wo + U)$ . 



(31) 



with the abbreviations 
A 



1 



W 2 - 



A 
2~ 



A = — 



9x = 9x 



Vx{l + tlx) 



A(l-Tb)(l + 3>fa) 
8 1 + 77, 



»7a 



2% 



(32a) 
(32b) 

(32c) 



It holds for the normal state: x = 1, a;' = 2, A = A, 
5 = <5 and 771 =772 = 1; for the blue superradiant state: 
a; = 1,2/ = 2, A = A, 6 = S and r/i = (gi/gi, c ) 2 ', and 
for the red superradiant state: x = 2,x' = 1, A = A — 5, 



5 = -6 and ry 2 = (WflVi)' 



IV. PHASE TRANSITIONS 

Comparing the ground-state energies of the states 
we found in the last section, we can derive the zero- 
temperature phase diagram. As mentioned before, the 
normal state is only stable for couplings g\ < g\ tC and its 
energy is independent of both coupling strengths g\ and 
g 2 . We also observed that the energy of the blue super- 
radiant state is always less than the energy of the normal 
state. However, the blue superradiant state is stable for 
.9i > 9i.c only. In addition, by comparing the energies 
of the blue (24 1 and the red (27 1 superradiant state, we 



see that only for g 2 > g 2 , C2 the red superradiant state is 
stable. Furthermore, in this parameter regime its energy 
is always smaller than the energy of the normal state (cf. 
Eq. pfk with g 2 = g2,c 2 )- 



A. Phase diagram 

From this discussion we derive the phase diagram 
which is shown in Fig. [3] It consists of three phases: 
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Figure 3: (Color online) Phase diagram for A > 5 > show- 
ing the three different phases: the normal and the blue and 
red superradiant phase. The (symmetric) normal phase is 
defined by $2 = ^3 = fi = <P2 — 0. In the (symmetry- 
broken) blue superradiant phase ^2 = 0, $3 / 0, 1^1 / and 
952 = 0. Finally, in the (symmetry-broken) red superradiant 
phase ^2 7^ 0, ^3 7^ 0, ipi — and ip2 7^ holds. The phase 
transition from the normal to the blue superradiant phase is of 
second order (red dashed line), whereas the phase transition 
from the normal to the red superradiant phase and between 
the two superradiant phases is of first order (red solid line). 
The normal state is meta stable in the region of the red su- 
perradiant phase as long as gi < gi. c . 



one normal phase for couplings g\ < gi iC and g 2 < 32,021 
one blue superradiant phase for couplings gi > <? l c and 
.92 < ff2,c(<?i)j and lastly one red superradiant phase for 
couplings g\ < g\ tC (g 2 ) and g 2 > g 2 , C2 - If both cou- 
plings are at criticality, gi = g± tC and g 2 = .92, C2 j all three 
phases coexist, i.e. there is triple point in the phase dia- 
gram. Here, 51,0(52) and g 2 , c (gi) parameterize the same 
curve, which represents the phase boundary between the 
two superradiant phases (see Fig. ptf. Both 51^(32) and 
<72,c(<7i) are given by the condition that the energies of the 
blue p4| and the red (p7|) superradiant state intersect, 
i.e 



both can be obtained by setting the Eqs. (24) and 



(27 1 equal. For gi^ c (g 2 ) we obtain after several algebraic 



transformations 



— 2 / \ 2 1 w l 

9iA92) = g 2 



2 UJ 2 



g2, Cl 
92 




(33) 



In the limit 6 — > 0, the phase bounda ry flattens to a 
straight line, lim^ §1,0(92) = 

The order of a phase transition is defined by the non- 
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Figure 4: (Color online) Ground state energy h/°\ ground- 
state occupation ip\ of the first (n = 1) and the second (n = 
2) bosonic mode and the ground-state occupation vE^ of the 
single particle energy levels (n = 1,2,3). Numerical values: 
A = Wi = 1, 8 = 0.75, (^2 = 0.25 (on resonance). 



analytic behavior of a thermodynamic potential |19| . 
In the case of zero temperature, the ground-state en- 
ergy represents a thermodynamic potential and hence its 
derivatives give the order of the phase transition. The 
ground-state energy of the normal state is E\ irrespective 
of the couplings 51 and 32. Hence, all derivatives with re- 
spect to <7i and 52 vanish. Comparing this result with the 
first and second derivatives of the ground-state energy of 
the blue (24) and the red (27 1 superradiant state, we see 



that the phase transition from the normal phase to the 
blue/red superradiant phase is of second/first order. The 
ground-state energy is shown in Fig. [4] 

In addition, the parameters ty r (r = 2,3) and tp n 
(n = 1, 2) also give evidence for the phase transition and 
can be interpreted as order parameters. An order param- 
eter is continuous for second-order phase transitions and 
discontinuous for first order phase transitions [19|. This 
behavior is visible in Fig. [4] The order parameters are 
zero in the symmetric (normal) phase and are finite in 
the symmetry-broken (superradiant) phase. The corre- 
sponding symmetry is the parity symmetry (see Sec. |TT|. 
In the blue/red superradiant phase the parity symmetry 
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Figure 5: (Color online) Excitation energies e„ jCT from 
Eqs. |30l, pll. Numerical values: A = wi = 1, 5 = 0.75, 



LU2 — 0.25 (on resonance). 



corresponding to the parity operator II1/II2 (see Eq. (|3|) 
is broken, since e.g. in the blue superradiant phase for 



finite ipi the operator c\ c\ in the Hamiltonian (18 1 is 
not invariant under the symmetry transformation Hi: 
Hi c\ ci fl\ = c\ ci + VJ7 (pi(c{ + ci) +Af<pf. 

Both, the phase transition and the order of the phase 
transition can also be deduced from the excitation ener- 



gies. The excitation energies from the Eqs. (30) and (31 1 
are shown in Fig. [5] At the phase transition at least 
one of the excitation energies either tends to zero or is 
discontinuous. The first case corresponds to a second- 
order, the latter case to a first-order phase transition. 
The second-order phase transition can be read off the 
excitation energy e\— which is zero for g\ = g\ tC and 

92 < 92,c 2 - 



Finally, we note that the phase transition from the nor- 
mal to the blue superradiant phase is in accordance with 
the superradiant phase transition in the Dicke model |15| , 
i.e. it is of second order and one (atomic) branch of 
the excitation energies tends to zero at the phase transi- 
tion. The discontinuity of the order parameters and the 
first derivative of the ground-state energy at the phase 
transition between the normal and the red superradiant 
phase scales with y/6. Thus, this first-order phase tran- 
sition becomes continuous in the limit S — > 0. However, 
the phase boundary between the two superradiant phases 
persists to be a first-order phase transition in this degen- 
erate limit. This is also the case in the limit of large 
couplings, 51/31, c, 52/52, c 2 00. 
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B. Dark state 

Due to the interaction of a quantum system with its 
environment decay processes within the quantum system 
occur. Eigenstates of the Hamiltonian which are unaf- 
fected by these decay processes are called dark states. In 
our model, a dark state is a many-body state which docs 
not radiate, i.e. a state where the occupation in either 
of the bosonic modes is zero. This condition is satisfied 
if the two parameters ip n are zero. From the Eqs. (20) 



we see, that the normal state, with \&2 = ^3 = 0, is a 
trivial dark state. In general, it suffices to set ^3 = 
for a dark state. Applying this condition to the Eq. (21 1, 
we can identify a dark state for 5 = only, i.e. for two 
energetically degenerate ground states. In the thermody- 
namic limit the coherence of this dark state is given by 
(A-^) = Afipi X I , 2) and is therefore finite apart from the 
two trivial cases -ipi = or $2 — 0. The energy of the 



dark state is simply h 



(o) 

Dark 



We obtain the excitation energies for the dark state 



' (2) 

by diagonalizing h y m ' =1 from Eq. (18 I. For any given \& 2 , 
these energies can be computed from the characteristic 
equation 



2ffi VW^i A 2g 2 ^V^A' 



= 0, 







w, — e 



det ( 2 i?i ipi^/ui A 
K 2g 2 *2V^2 A 

(34) 

where det M is the determinant of the matrix M. The 
characteristic equation is readily solved in the case of 
two-photon resonance (u)x — 102 = A = w), yielding the 
energies 



£0 = 0, Ei= 0J, 



£2 ,± = \^ ±2uJ g m - * 



2 I 



(35) 
(36) 



where the additional zeroth mode £0 stems from the limit 

The parameter vf^ is arbitrary and can range from zero 
to one. For a given we find by analysis of the Hessian 
matrix of ^Li> that this dark state is meta-stable if the 
inequality 



(ffl/ffl,c) 2 (l-*l) <l-(ff2/52,c) 2 ^ 



(37) 



is satisfied. Otherwise this dark state solution is unsta- 



ble. In Eq. (37) we have introduced the critical coupling 
strength g% c = y/A w 2 /2. 

We emphasize that the dark state exists for 6 = 



only. By inspection of the inequality ( 37 1 , we make the 



following statements: First, the dark state is stable for 
9i < 9i,c or <?2 < ff2,c only. Furthermore, if both coupling 
strengths fulfill g n < g n ,a i- e - m the normal phase, both 
and ^>2 can range from zero to one. On the other 



hand, if 52 > 92, c an d .91 < gi, c , then ^2 is restricted to 
the interval [0, ^2, max], where ^2, max > is given by the 
inequality (37 1. Correspondingly ipi is restricted to the 

interval [?/>i, m in, 1], with Vl.min given by yjl- ^l.max- 
An analogue argument can be given for the case g\ > g\ >c 
and g2 < 32, c, where ipi and \&2 are interchanged. 

For couplings g\ ^> g\ >c and g2 < g2.c, inequality (37 1 
restricts the order parameters to w and ^2 ~ 1, i-e. 
only the second single-particle level is macroscopically 
occupied. On the other hand, for couplings g 2 3> g2, c 
and gi < g\^ c only the first single-particle level is macro- 
scopically occupied, i.e. i\>y k> 1 and $2 ~ 0. This 
'counterintuitive' behavior is reminiscent of the STIRAP 
scheme [TJ]. In contrast to the STIRAP scheme, the 
actual values of the populations ipi and ^ 2 in this dark 
state are not defined by the coupling strengths g± and 172 , 
but rather by the preparation of the system. Thus, the 
system cannot be driven coherently from a state with all 
particles occupying the first single-particle energy level 
|1) to a state where all particles occupy the second single- 
particle energy level |2) just by changing the couplings. 

In addition, we note that in this dark state the 
mode Eq = in direction of ^2 of the energy surface 
ft( )($ 2 ,xp 3 ) is trivially massless (cf. Eq. d35j). There- 
fore, tiny fluctuations can easily excite this dark state 
along the direction of ^» 2 , making the state eventually 
unstable. This instability is visualized in Fig. [6] 



V. CONCLUSIONS 

We have analyzed an extension of the well known Dicke 
model from two to three-level particles. By means of 
a Holstein-Primakoff transformation we have identified 
three stable states in the thermodynamic limit: a nor- 
mal, a blue superradiant, a red superradiant state. At 
zero temperature these states correspond to three ther- 
modynamic phases, which we have arranged in a phase 
diagram. The phase transition between the normal and 
the blue superradiant phase is of second order and all 
other phase transitions are of first order. We have also 
shown that a state with both superradiant states coexist- 
ing is not stable. A dark state with zero occupancy of the 
third single-particle level exists for 6 = only. However, 
this dark state is not stable. 

As in the original Dicke model, the same experimental 
difficulties arise in our extended Dicke model, i.e. reach- 
ing the critical coupling strength is challenging as well. 
Hence, using three level atoms has no advantage over the 
use of two level atoms. 

However, we expect that similarly to the Dicke model 
and its realization in the experiments of Baumann et al. 
[TT] there should be experimental manageable systems, 
which can theoretically be described by an effective 
Hamiltonian of the form presented here. In the case 
of the experiments in Ref. |llj , this might be achieved 
by coupling a Bose-Einstein condensate to an additional 
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cavity mode. Furthermore, an even richer phase diagram 
with additional superradiant phases could be generated 
in such a system. 

Considering a cold quantum gas in an optical lattice, 
a characteristic feature of our extended Dicke model es- 
pecially in the degenerate limit, 5 0, could appear. 
In this regard, we have an extension of a system pro- 
posed by Silver et al. |20| in mind. There, it was shown 
that a two-band zero-hopping Bose-Hubbard model cou- 
pled to a cavity light field can be written as an effective 
Dicke model. If one superposes a superlattice of twice 
the wavelength of the original lattice, and couples the 
superlattice to two independent cavity light fields, this 
extended Bose-Hubbard model can be mapped to our 
extended Dicke model with 5 = 0. Since in experiment 
one has an extensive control over the parameters of cold 
quantum gases, the observation of superradiant phases 
should be feasible. 



l-i l-i 



Figure 6: (Color online) Phase diagram (a) as in Fig. [3] and 
energy surfaces /r ' (b-d) in the degenerate (5 — > 0) limit. 
In (a), the solid and the broken red lines denote first and 
second-order phase transitions, respectively. The red lines in 
(b-d) visualize the inequality (37 1, i.e. the line where the 
dark state with $3 = is stable. However, as the red line 
is flat, fluctuations along the ^-direction transform stable 
states on the red line into unstable states outside the red 
line. Eventually, these states decay to superradiant states 
with * 3 ^ 0. 
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